Star formation rate density and the stochastic background of gravitational waves 
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There is in the literature a number of papers addressing the stochastic background of gravitational 
waves (GWs) generated by an ensemble of astrophysical sources. The main ingredient in such studies 
is the so called star formation rate density (SFRD), which gives the density of stars formed per unit 
time. Some authors argue, however, that there is, in the equation that determines the amplitude 
of the stochastic background of GWs, an additional (1 + z) term dividing the SFRD, which would 
account for the effect of cosmic expansion onto the time variable. We argue here that the inclusion 
of this additional term is wrong. In order to clarify where the inclusion of the (1 + z) term is really 
necessary, we briefly discuss the calculation of event rates in the study of GRBs (gamma ray bursts) 
from cosmological origin. 
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I. STOCHASTIC BACKGROUND OF 
GRAVITATIONAL WAVES 

Because of the fact the gravitational waves (GWs) are 
produced by a large variety of astrophysical sources and 
cosmological phenomena it is quite probable that the 
Universe is pervaded by a background of such waves. 

We have recently shown jl| that the dimensionless 
amplitude of the stochastic background of GWs produced 
by an ensemble of sources is given by 



h 2 



1 

t'obs 



B dR, 



(1) 



where /i S mgic is the dimensionless amplitude produced by 
an event that generates a signal with observed frequency 
fobs! and dR is the differential rate of production of GWs. 

Eq.lpJ is in fact a shortcut to the calculation 
of stochastic background of GWs. An interesting 
characteristic of this equation is that it is not necessary 
to know in detail the energy flux of the GWs produced 
at each frequency. If the characteristic values for the 
dimensionless amplitude and frequency are known and 
the event rate is given it is possible to calculate the 
stochastic background of GWs produced by an ensemble 
of sources. 

In particular, for the case of a background produced 
by an ensemble of black holes, the differential rate reads 



dR^n = p Sf (z)——(j)(m)dmdz, 
dz 



(2) 



where p» [z] is the star formation rate density [hereafter 
SFRD (we refer the reader to the Appendix, where it 
is shown how the SFRD is obtained from observations); 



in M0 yr~ x Mpc ], dV is the comoving volume element 
and 4>( m ) is the initial mass function (IMF). In a few 
words, the 4>(m) dm represents the number of stars per 
unit mass in the interval [m, m + dm\. The normalization 
of the IMF is obtained through the relation 



m (fi(m) dm = 1, 



(3) 



and 
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where it is usually considered that mi = O.1M 
m u = 125M (see, e.g., Ref. Q). 

Some authors argue that there is an additional (1 + z) 
term in the equation for the differential rate dividing 
p*(z), which would take into account the effect of cosmic 
expansion onto the time variable. 

In fact, there are two different points of view in the 
literature concerning the inclusion or not of the factor 
(1 + z). The first group does not include this factor 
dividing the SFRD (see, e.g., Refs.0,HI3). 0n thc othcr 
hand, the second group argues that it is necessary to 
include the factor (1 + z) to account for the time dilation 
of the observed rate by cosmic expansion, converting a 
source-count equation to an event-rate equation (see, e.g., 
Refs.0, 0, |(| 0). We show here why such an additional 
term does not exist. 

In order to address properly this issue, we should note 
firstly that the determination of the amplitude of the 
GWs is obtained from the flux received by a detector 
(.Few)- In this way, one divides the luminosity of the 
sources by And^, where c?l is the so called luminosity 
distance. The luminosity obviously refers to the source 
frame. Since the luminosity depends on the differential 
rate, given, for example, by Eq.Q, the latter also refers 
to the source frame, therefore it is not necessary to 
redshift it. It is worth stressing that any necessary 
redshifting is taken into account in thc definition of d\^. 

Another way to be convinced that the (1 + z) does 
not enter into the calculation of the differential rate of 
production of GWs, cLRbh, is to derive Eq.QJ through a 
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procedure completely different from that used in Ref.pJ. 

Let us write the specific flux received in GWs at the 
present epoch as (see, in particular, Eq. (15) in Ref. 8] 
and section 12.1 in Ref.0) 



Id dv 
Awd? dv ohs 



dV, 



where 



dL v 



(4) 



(5) 



is the comoving specific luminosity density (given, e.g, 
in ergs -1 Hz -1 Mpc -3 ), which obviously refers to the 
source frame. 

As discussed in Refs.jl Q, the above equations are 
valid to estimate a stochastic background radiation 
received on Earth independent of its origin. In the 
present paper l v can be written as follows 



'" = ^ ^dv^ P * ( z W TO )^ TO ' 



(0) 



where dEcw/dv is the specific energy of the source. Note 
that in the above equation p*(z) refers to the source 
frame, therefore, there is not the putative (1 + z) factor 
responsible to the time dilation. 

Thus, the flux F v {y \^) received on Earth reads 



F v {v oha ) 



1 dE, 



GW 



dv 



And? dv dv bs 



p*(z)4>(m)dmdV. (7) 



Using Eq.J2J it follows that 

1 dEcw dv 



And? dv dvoXjt 



■dRsR- 



(8) 



Note that in the above equation, what multiplies cLRbh is 
nothing but the specific energy flux per unity frequency 
(in, e.g., erg cm -2 Hz -1 ), i.e., 



1 dE, 



GW 



dv 



^ obs '' 4irdf dv dv hs 



(9) 



(see, e.g., Ref.[To|'). 

On the other hand, the specific energy flux per unit 
frequency for GWs is given by Ref.[TT| 



Jv{Vobs) = ~2Q BK ' 



(10) 



Also, the spectral energy density, the flux of GWs, 
received on Earth, F v , in erg cm -2 s -1 Hz -1 can be 
written from Refs. [l%ll3| as 



F v (v \,^) — "^T^BG^obs- 

From the above equations one obtains 

h?3G = / ^BH^BH, 

^obs J 



(11) 



which is nothing but the Eq. QJ. 

One could argue that what we derive here is 
nothing but a rehash of a well known equation for 
the calculation of electromagnetic backgrounds, found 
in textbooks on cosmology, applied to gravitational 
radiation backgrounds. But, since some authors 
in the GW community are erroneously introducing 
extra redshift factors in their calculations, it is worth 
presenting clearly how to calculate GW backgrounds of 
cosmological origin. 

Also, these authors do not base their calculations 
either on the formulation of electromagnetic backgrounds 
presented in textbooks on cosmology or on an equation 
such as the one derived here. They inappropriately base 
their calculations on event rate equations, in which the 
time dilation needs to be considered. As a result, they 
consider that in other calculations involving the SFRD 
the redshifting should be taken into account. 

It is worth stressing that the equation for the energy 
flux of GWs, we use in our derivation, is not obtained by 
adapting the equation for the Poynting flux by replacing 
the electric field with the GW amplitude, with the 
corresponding factors of G and c. This equation comes 
from general relativity. 

One could argue when it is necessary to take into 
account the redshifting related to the time dilation. We 
now discuss a case where we must include the factor 
(1 + z) to take into account the time dilation of the 
observed rate. 

The case is related to the Gamma Ray Bursts 
(hereafter GRBs), which are short and intense pulses 
of 7— rays, which last from a fraction of a second to 
several hundred seconds. It is worth mentioning that 
there are in the literature many papers dealing with the 
GRB event rates of cosmological origin, namely, Refs. 
El El El 113 > among others. 

The event rate for the GRBs of cosmological origin is 
related to the SFRD. The main reason to include the 
(1 + z) factor in this case is that we observe GRB events 
over a fixed time window Ai b s , which corresponds to 
Atobs/ (1+- 2 ) m the source frame. Thus, the total number 
of GRBs can be written as (see, in particular, Ref-E3) 



N{> z) 



1 , At ohs dV , 

^GBB{z) 7TT ^—dz , 



(1 



(12) 



where dV/dz is the comoving volume element and 
V'grb(z) is the number of GRB events per comoving 
volume per unit time, which is proportional to p\{z). 

Generally speaking, in any event rate of cosmological 
sources, which involves a formation rate as a function of 
redshift or the like, the time dilation must be applied. 

We see that the definition presented in Ea. l|12[l is 
completely different from the definition presented in 
Eq.Q). The misuse of the time dilation of the SFRD in 
the calculation of the energy flux is due to the fact that 
this involves a time rate, as the event rate of cosmological 
sources does. But, one has to bear in mind that the 
luminosity distance, that links the luminosity (source 
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frame) and the energy flux (observer frame), is defined 
in such a way that any redshifting, including the time 
dilation, is implicity in its definition. 

In the next section, we reinforce the argumentation 
against the inclusion of the extra (1 + z) term recalling 
how luminosity, flux and the luminosity distance are 
related to each other. 



II. LUMINOSITY, FLUX AND LUMINOSITY 
DISTANCE 

Although the material present in this section is 
standard, appearing in any textbook of cosmology, it is 
worth having a look at the relation between luminosity, 
flux and luminosity distance in order to recall that any 
redshifting, in particular the time dilation, is already 
taken into account. 

An isotropic emission of photons (or gravitons in the 
present case) emitted by a source pass through spheres 
surrounding the source and, in the absence of expansion, 
the detected flux is exactly equal to the fraction of the 
area of the sphere surrounding the source and covered by 
the detector, that is, dA/Aird^ (where dA is the area of 
the detector) times the luminosity L of the source. In 
this case, g£l is simply the distance to the source (at this 
point we refer the reader to some textbook on cosmology 
as, e.g., iEEl). _ 

When the expansion of the Universe is taken into 
account, this consideration is slightly modified. In this 
case, we can use the Friedmann - Robertson - Walker 
(FRW) metric as being centered on the source; because of 
homogeneity, the comoving distance between the source 
and the observer is the same as we would calculate when 
we place the origin at the observer's location. 

Setting t = constant and r — constant in the FRW line 
element, then we obtain 



ds 2 = -r 2 a 2 (d9 2 



sin 



(13) 



where this equation represents the line element on the 
surface of a Euclidean sphere of radius ar, and a = a(t) 
is the scale factor of the Universe. 

The total energy emitted by a source per unit time, at 
the epoch t\, over the bandwidth (v, v+ A v) is given by 



dL = L J(v) du, 



(14) 



where J(v) is the intensity function. 

In the case of a source emitting isotropically, when its 
light reaches us it is distributed uniformly across a sphere 
of coordinates r = r\ centered on the source. However, 
we should note that for the source situated at comoving 
coordinate r — r\ the photons (or gravitons) emitted at 
time t\ will be detected at comoving coordinate r = at 
time to. 

Note that because of the redshift, the arriving photons 
(or gravitons) with frequencies in the range (v bsi v bs + 



A f bs) left the source in the frequency range [f bs(l + 
z), (v oha + Av ohs )(l + z)}. 

Thus, the quantity of energy that leaves the source (at 
r = ri ) between the times t\ and t\ + Ati in the above 
frequency range is simply obtained by 



L J(v ohs [1 + z]) Ai> ohs (l + z) Ati. 



(15) 



In order to find the number of photons (or gravitons) that 
carry the above quantity of energy, we should divide Eq. 
(|15fl by the energy (1 + z) /i^obs, where h is the Planck 
constant. That is 

LJ(v ohs [1 + z]) Aii . Mfi . 

5N = Av ohs . (16) 

hv h s 

At the time of reception, these photons (or gravitons) are 
distributed across a surface area A-kt\ a 2 (to) and they are 
received over a time interval (to, t Q + Ai )- Thus, the 
number of particles received by a detector per unit area 
held normal to the line of sight and per unit time is given 
by 



LJ(v ohs [l + z\) Aii 



hv ohs A t 47rr 2 a 2 (i ) 



A^obs- 



(17) 



When received, each photon (or graviton) has been 
degraded in energy by the factor 1/(1 + z). Thus, the 
photons (gravitons) now has the energy hi> b s . If we 
multiply Eq. i|17[l by this factor, we obtain the quantity 
of photons (gravitons) received per unit time by us across 
unit proper area held perpendicular to the line of sight 
to the source, and over a bandwidth (^ bs, ^obs + A fobs)- 
Denoting this quantity F^(v \, s ) A f b s , we have 



F v (v ohs ) A i/ obs = LJ(v ohs [1+z]) 



1 



Ati 

Ai^ 4irrl a 2 (to) 



Au 



obs ■ 



(18) 

We should note that in the above equation Ati/At = 
1/(1 + z). This is just the (1 + z) factor claimed 
by some authors that should divide the SFRD to take 
into account the time dilation of the observed rate by 
cosmic expansion. As seen below this factor is already 
incorporated in the definition of the luminosity distance. 
Thus, finally we obtain 



Fv(v bs) 



LJ(u ohs [1 + z]) 
4nr 2 a 2 (t )(l + z)' 



(19) 



where Fi/^obs) is the flux density. 

In a few words, the above result can be understood 
as follows: The photons (gravitons) from the source 
pass through a sphere of proper area Airr 2 on which 
the observer sits, where r\ is the comoving distance. 
Moreover, redshift affects the flux density in four ways: 
photon (graviton) energy and arrival rates are redshiftcd, 
reducing the flux by a factor (1 + z) 2 ; opposing this, 
the bandwidth A v is reduced by a factor (1 + z), so 
the energy flux per unit band goes down by one power 
of (1 + z); finally, the observed photons (gravitons) at 
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frequency ^ bs were emitted at frequency v Q hs (1 + z). 
Thus, the flux received on Earth from a cosmological 
source is the luminosity of the source divided by the total 
area and divided by (1 + z) as given in Eq. i|19|) . 

The bolometric flux can be obtained from the 
integration of Eq. (|19J) over all frequencies. The result is 



^ b0l ~ 4^«2(* )(l + z)2' ^ 

where the quantity c2l = r\ <x(io) (1 + z) is the so called 
luminosity distance. It is worth stressing that Eq. H20|) 
is a direct consequence of the conservation of energy. 



community, or any one who calculates cosmological 
event rates, should keep using the time dilation in their 
calculations. 
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III. CONCLUDING REMARKS 

The (1 + z) dividing the SFRD, which some authors 
claim would account for the time dilation of the observed 
rate by cosmic expansion, converting a source-count 
equation to an event-rate equation (see, e,g, 0,IEIEI3l)j 
is nothing but the arrival rates, which is implicit in the 
luminosity distance. In particular, Refs. 2, 3] implicitly 
corroborate with the argumentation given here. 

It is important to stress that for the cosmological 
background of GW models we studied in Ref.j^, the 
inclusion of the (1 + z) term would reduce by a factor 
of ~ 3 — 5 the signal-to-noise ratio (SNR) predicted for 
the LIGO observatories. For some models, this would 
mean to predict a non detection of such a stochastic 
background of GWs. 

We hope this brief report contributes to clarify the 
reader why in the calculation of the amplitude of 
a stochastic background of GWs, from cosmological 
sources, it is not necessary to redshift the SFRD as 
some authors are arguing for; and also why the GRB 



APPENDIX: THE STAR FORMATION RATE 
DENSITY 

The SFRD is inferred from observations of the light 
emitted by stars at various wavelengths. In particular, 
|2l| investigated the galaxy luminosity density of rest- 
frame ultraviolet (UV) radiation up to z ~ 4, and they 
converted it into the SFRD. The rest-frame UV light is 
considered to be a direct tracer of star formation because 
it is mainly radiated by short-lived massive stars. 

These observable samples are flux-limited, and thus 
the intrinsic luminosity of the faintest objects in the 
sample changes with redshift. This incompleteness of 
the samples is corrected by using a functional (Schechter 
function) to the luminosity function obtained from the 
observations themselves. Then, the conversion from 
luminosity density to SFRD generally relies on stellar 
population models, an assumed star formation history 
and a specific choice for the IMF (generally, it is 
considered the Salpeter IMF). 

In this way, the SFRD is mainly derived from the 
observed luminosity at high redshift, that is, p*(z) oc L. 
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